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ABSTRACT 



An analytical investigation of laminar free convection 
from a semi - inf inite vertical plate is undertaken with the 
following special conditions imposed. The plate is main- 
tained at a constant temperature greater than that of the 
ambient fluid. At an arbitrary distance above the leading 
edge a constant rate of blowing or suction is initiated 
extending uniformly along the remainder of the plate. The 
solution was obtained using the method of matched asymptotic 
expansions whereby an inner series described the velocity 
and temperature profiles near the discontinuity in the 
vicinity of surface and an outer series approximated the 
behavior at greater distances from the plate. The inner 
series was also used to determine the heat transfer charac- 
teristics at the surface of the plate. The ambient fluid 
was assumed to be air with a Prandtl number of 0.72. 
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I. INTRODUCTION 



The study of free convection has been the subject of 
ongoing research over the past century both from an analyti- 
cal and experimental point of view. An excellent overview 
of the progress made in this area can be found in Ede [1] . 

Two outstanding expositions of the fundamentals of heat 
transfer by the mechanism of free convection are presented 
by Gebhart [2] and Moore [3]. 

The technique of solving boundary layer problems by 
employing similarity analysis is discussed in general by 
Hansen [4] and specific applications of the method are pre- 
sented by Yang [5]. 

The problem discussed in this paper is that of an analyti- 
cal investigation of the effects of blowing and suction 
applied at a specified distance above the leading edge of a 
uniformly heated semi - inf inite vertical plate on the heat 
transfer characteristics at the surface of the plate. 

An early study of the effects of mass transfer on free 
convection was conducted by Eichhorn [6] who showed that 
similar solutions are possible for blowing or suction rate 
distributions which vary as the distance from the leading 
edge of the plate raised to a power related to the variation 
in surface temperature. 
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Since the prescribed surface temperature itself varied as 
a power of the distance from the leading edge, this proved to 
be a somewhat restricted case of limited physical application. 
The case of a uniformly heated plate with constant blowing or 
suction along the entire extent of the surface was succinctly 
treated by Sparrow and Cess [7] who extended their analysis 
to include independent power law variations of both surface 
temperature and mass transfer. Their analysis which was based 
on an asymptotic series expansion of the stream function and 
temperature distribution evaluated the first order effects of 
blowing and suction on free convection. An extension of this 
work to include possible variations in the fluid properties 
was conducted by Parikh, Kays and Bershader [8] who employed 
a numerical finite difference scheme to solve the resulting 
equations. Their numerical predictions were in close agree- 
ment with interferometric data obtained through their own 
experiments. The results of these authors and especially 
those reported by Sparrow and Cess have been of great value 
in affording a comparison with results obtained in the 
present study. A comprehensive investigation of the effects 
of suction and blowing on heat transfer and wall skin fric- 
tion which included regions on the plate at large distances 
from the leading edge was conducted by Merkin [9] . In that 
analysis a step-by-step numerical solution of the boundary 
layer equations begins at the leading edge and proceeds 
systematically up the plate until the asymptotic velocity 
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and temeperature profiles are obtained to whatever accuracy 
is required. 

There have been several additional studies in the area 
of free convection past a vertical plate with various condi- 
tions of wall temperature or heat flux and imposed variations 
in transpiration at the surface of the plate. An analysis 
which transformed the laminar boundary layer equations re- 
sulting from arbitrary distributions in both wall temperature 
and blowing into ordinary differential equations with variable 
coefficients and constant boundary conditions was presented 
by Vedhanayagan, Altenkirch, and Eichhorn [10]. In the area 
of unsteady free convection where either the plate surface 
temperature or transpiration rate varies with time there have 
been at least four notable studies during the past twenty 
years. Nanda and Sharma [11] undertook analytical inves- 
tigations and determined that similarity solutions were 
possible for the following two cases: surface temperature 

varying as some power of time and suction velocity varying 
as time to the minus one-half power; and the second case 
where surface temperature is an exponential or periodic 
function of time and the suction velocity is constant. The 
effects of mass diffusion on the unsteady free convective 
flow past a simi- infinite porous plate with constant suction 
were studied through mathematical analysis by Soundalgekar 
and Wavre [12] . With the assumption that the plate temperature 
oscillated in time about a constant mean, approximate 
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solutions were derived for the mean flow, transient flow, the 
amplitude and phase of the skin friction and the rate of heat 
transfer. The effects of the nondemens ional parameters Grashof 
number, Prandtl number, Eckert number and Schmidt number were 
also reported. The most general unsteady situation where both 
wall temperature and suction velocity are periodic functions 
of time was considered by Lai [13] . All of the previous 
studies cited have neglected the effects of viscous dissipation; 
however, this effect is not insignificant when the natural 
convection flow field is very large or the flow occurs in the 
presence of extremely low temperatures or very high gravity 
fields. The combined effect of viscous dissipation with the 
suction velocity and plate temperature oscillating at the 
same frequency on free convection past a vertical plate was 
analyzed by Soundalgekar and Pop [14]. In their study expres- 
sions for the mean velocity, transient velocity, transient 
temperature profiles, amplitude and phase of skin friction 
and the rate of heat transfer were derived and the results 
compared to those reported for the case of constant suction 
velocity. 

A review of the literature in the field of free convection 
on a vertical plate clearly indicates a desire on the part 
of researchers to investigate perturbations of the thermal 
and momentum boundary layers resulting from the introduction 
of temperature or velocity discontinuities along the surface 
of the plate and their ultimate effect on skin friction and 



16 



heat transfer at the wall. The specific nature of the dis- 
continuity will determine the manner in which the perturbation 
effects propagate through the boundary- layers and; therefore, 
will often indicate a preferred method of analytical solution. 
Although there are several possible mathematical techniques 
available to researchers in this field, two fundamentally 
different approaches predeminate, each with numerous variations. 
An example of the finite difference method of solution to the 
boundary- layer equations resulting from a discontinuity in the 
surface wall temperature can be found in Na [15] and Meena and 
Nath [16] . This technique has gained in acceptance with the 
advent of larger and faster digital computers. The other 
major technique involves the expansion of both the stream 
function and temperature profile in an asymptotic series. 

The original proponent of this method was Goldstein who pre- 
sented the details of the asymptotic series technique, 
commonly known as the Goldstein series, in [17]. After re- 
ducing the set of governing boundary- layer equations by the 
introduction of a stream function such that: 

u = 37 and v = -57 

two new variables are introduced which are of the following 
general form: 

1) In the streamwise direction: £ = x 11 

- 1 /n 

2) In the direction normal to the plate: n = yx 

A series is then formed in powers of £ for the stream function 
and temperature distribution with undetermined coefficient 
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functions f ^ respectively which depend only on n . The 
problem is then finally reduced to one of determining the 
coefficient functions; the solution of a fifth order set of 
linear coupled ordinary differential equations for each pair 
of coefficients subject to the appropriate boundary conditions. 
This technique is especially powerful when employed to inves- 
tigate conditions at or near the surface of the plate. In 
order to examine the laminar free convection wake above a 
heated vertical plate Yang [18] successfully employed this 
technique to obtain detailed velocity and temperature profiles 
in the immediate vicinity of the trailing edge and then 
expanded these results to the rest of the wake by an integral 
solution. Another example of this technique is found in 
Kelleher [19] where a step discontinuity in the surface 
temperature imposed at a finite distance above the leading 
edge of the plate caused a perturbation of both the thermal 
and momentum boundary- layers . An inner expansion of the 
Goldstein type was employed to obtain the behavior of heat 
transfer near the surface of the plate and an outer asymptotic 
series expansion was constructed to determine the velocity 
and temperature profiles at greater distances from the plate. 

A similar approach is employed to analytically investigate 
the present problem where a discontinuity in the velocity 
boundary condition is propagated through the momentum and 
thermal boundary- layers . Laminar flow is assumed as are 
constant transport properties of the ambient fluid. The 
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primary aim of the study is the determination of any en- 
hancement or degradation of the heat transfer rate at the 



surface of 
blowing or 



the plate compared to the same 
suction over the entire length 



plate without 
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II. ANALYSIS 



A. COORDINATE SYSTEM 

The problem considered in this paper is that of laminar 
free convective flow from a uniformly heated semi - inf inite 
vertical plate with leading edge oriented downward in a 
standard gravitational field. A constant rate of blowing or 
suction is applied at the surface beginning at an arbitrary 
distance, L, above the leading edge. A rectangular cartesian 
coordinate system with origin as indicated in Figure 1 was 
employed in the analysis. As is indicated throughout the 
paper various transformations of the basic coordinates are 
invoked to facilitate the solution of the problem. 

B. GOVERNING EQUATIONS 

The analysis is based on the standard Boussinesq free 
convection approximation which assumes that only density va- 
riations giving rise to buoyant forces are considered with 
the effects of viscous dissipation neglected. With these 
standard assumptions the two dimensional equations of con- 
tinuity, momentum, and energy boundary - layer form are: 

* If . o (l) 

3x 3y 

u + v ^ = v ^ + g0(T - TJ (2) 

3x 3y 3y 2 
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Figure 1: Basic Coordinate System 



- 3T - 3T 3 2 T 

u — + v — = a 



( 3 ) 



3x 3y Sy^ 

where is the constant ambient fluid temperature. These 
equations must also satisfy the following physical boundary 
conditions : 



At 


y = 0 : 


u = 0 








d 

n 


x>0 






O 

II 

l> 


- L<x<0 






T = T 

w 




as 


y -* oo • 


u = 0 





(3a) 



T = T 

oo 

The non- zero normal velocity boundary condition at x = 0 and 
y = 0 introduces a discontinuity into the boundary- layer equa- 
tions which has a primary direct effect on the formation of 
the momentum boundary- layer and a secondary influence on the 
development of the thermal boundary- layer . A positive value 
for will represent the blowing of fluid away from the wall 
and a negative value for v will represent the suction of fluid 
through the plate. 

In order to facilitate the analysis, the variables in 
equations (l)-(3) are non-dimensionalized according to the 
following transformations: 



x = r 



y 



= y 



(4a) 

(4b) 
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u 



uL 

V 



(4c) 



v = 



vL 

v 



( 4d) 




(4e) 



where L is the distance above the leading edge where suc- 
tion or blowing commences, and is the constant wall 
temperature. The equations are further simplified through 
the introduction of the conventional stream function ip and 
by stretching y, u, and v as follows: 



mV 
Jh IrJ- 

X 

II 


(5a) 


= 0(401)^ = |1, 


(5b) 


- v(64Gr)'* s - - || 


(5c) 



where the Grashof number, Gr , is evaluated with respect to 
L, the length scale associated with the discontinuity. It 
should be noted that this length scale is not the conven- 
tional one used in free convection problems where distance 
from the leading edge is normally used to evaluate Gr. Now 
after substituting equations (S) into equations (l)-(3), the 
continuity equation is eliminated and the momentum and energy 
equations reduce to the following: 



4 . dUj . 4 M . + a 

3y* 3x3y* 3x 3y* 2 3y* 3 



( 6 ) 
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d li li 

3y* 3x 



, dip 3 9 

4 37 37 ^ 



1 3 2 0 

Pr 3y*2 



( 7 ) 



The non- zero normal velocity component imposed at the surface 
of the plate perturbs the governing equations (6) and (7) 
precluding the existence of a similarity solution. There- 
fore, an approximate solution to these equations will be 
obtained by means of inner and outer asymptotic series ex- 
pansions for both ip and 9. The inner expansions will approxi- 
mate local behavior near the surface of the plate in the 
vicinity of the discontinuity, and the outer expansions will 
model velocity and temperature profiles at greater distances 
from the plate. 

C. BOUNDARY AND INITIAL CONDITIONS 

The nondimens ional form of the governing equations given 
in (6) and (7) must satisfy the original physical boundary 
conditions as specified in (3a) which reduce to the following 
conditions in terms of the stream function ip and temperature 
distribution 0: 

At y* = 0 : u* = = 0 

v* = = ±e* , x> 0 

v* = = 0 , -L<x<0 (8a) 

0 = 1 

As y* + u* = 0 
0=0 
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where e * is related to the physical blowing or suction rate, 



v , by the following: 
w 

- v L j 

e* = ~v~ ( 64Gr )"^ 



(8b) 



Since this problem involves a perturbation of both the 
momentum and thermal boundary- layers which began developing 
at the leading edge, the initial conditions at the origin of 
the discontinuity, i.e., x = 0, must match these free con- 
vection profiles. For a semi- infinite vertical plate main- 
tained at a constant temperature the following ordinary 
differential equations describe unperturbed free convection 
behavior : 





F T 1 ? H 


H 3FF" - 


2F ' 2 + H = 


= 0 




i- H" 
Pr 


+ 3FH ' = 


0 




with 


the following 


boundary 


conditions : 


At 


J 

II 

O 


F ( 0) = 


F ' (0) ; H (0) = 1 


As 


n+°° : 


F'(co) = 


0; H(oo) = 


= 0 



(9a) 

(9b) 



( 10 ) 



where F is the similarity stream function and H is the 
dimensionless temperature distribution and primes indicate 
derivatives with respect to the free convection similarity 
variable defined as follows: 



n f 



= y 

(x+L) 



r Gr* J* 
l 4 J 



(ID 
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where Gr* is evaluated at x+L according to the free convection 
length scale which has its origin at the leading edge. It 
should be noted that all references or comparisons to the 
unperturbed free convection problem employ this translated 
length scale, whereas the length scale for the perturbed 
problem is based on L with origin at x = 0. At x =0 equation 
(11) reduces to ri£ = y* . The solution to the equations (9) 
are well known; however, an excellent summary of solutions 
for this and other more general cases can be found in Yang 



Since the initial length of the plate is unperturbed, the 
velocity and temperature profiles at x = 0 will be given by 
the following: 



Therefore, the initial velocity and temperature profiles 
originating at x = 0 must match those obtained from equations 
(12) . The details of expressing F and H in a form which will 
permit this matching condition to be mathematically formulated 
are presented in Appendix A. 



[5]. 



u = (4Gr)^ F' (y*) 
0 = H(y*) 



(12b) 



(12a) 
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III. METHOD OF SOLUTION 



A. REDUCTION OF THE GOVERNING EQUATIONS 

The nondimens ional forms of the governing equations for 
the problem under consideration are: 



3y* 


d 2 ip 

3x3y* 


- 4 


3^ 3 z ip 
3x 3y*2 


= 3y*3 + 9 


(6) 


dy* 


M - 4 
3x 


d^p 

3x 


36 

3y* 


i 3 2 e 

Pr ay*^ 


(7) 



In order to apply the Goldstein series method of solution to 
these equations suitable expansions of ip and 9 must be con- 



structed of the following form: 

1 

Y - e*U/<Di) a + <M c U k f k (n)] 
and 9 = 1 + <j\£ d U k g k ( n) ] 



(13) 

(14) 



where E, =<(> x, n = <J> 2 y*x, the <p' s are possibly functions of 
e*; and the a, b, c, d are constants to be determined. The 
following conditions must be satisfied which impose certain 
restrictions on the <j>'s and exponent constants: 

1) The u velocities from the free convection region 

(x<0) and from the perturbed stream function, i.e., 

u * = -Ptf, must match at x = 0 . 
dy 

2) The temperature profiles must match at x = 0. 

3) The highest order derivative of the coefficient func- 
tions must appear in the zeroth order expansions which 
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result upon substitution of equations (13) and (14) 
into equations (6) and (7) . 

4) The exponent constants in the expansions referred to 
in 3) above must be integers. 

The details of applying the conditions outlined above to 
determine the unknown parameters are presented in Appendix B 
along with the boundary conditions which must be satisfied 
by the coefficient functions f^ and g^. 

Using the results of Appendix B equations (13) and (14) 
are transformed into the following: 

ip = 4*£ 3 + ? 2 U k f k ] (15a) 

and 9 = 1 + £[£ k g k ] (15b) 

where £ = (jx) 1//3 and n = (jx) " 1 / 5 y* = y*/£ . (16) 

The momentum and energy equations (6) and (7) are transformed 
into : 

[(k + l)£ k f’] [£ k f£] - [(k + 2)£ k f k ] [£ k f-] 

- 4c*£[£ k f£] = [£ k f k 1 ] + £ + £ 2 [£ k g k ] (17) 

and [(k + l)£ k g k ] [£ k f£] - [(k + 2) £ k f k ] [£ k g£] 

- 4e * £ [ £ k g k ) = -p7[5 k g k ] US) 

where the repeated index k indicates the standard Einstein 
summation. By the process of equating coefficients of like 
powers of £ the following ordinary differential equations 
involving the coefficient functions are obtained. 
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Momentum Equation: 
Oth order: 



f! ' ' + 2£ f" - f ' 2 = 0 
0 0 0 0 



order : 


f T 


T T H 


H 2f 


£ f! - 


3f T f T 


+ 3 f 


M f 


= - 4e * 


£ n" 


- i 


(19) 




1 




9 


i 


0 i 




0 1 




0 




order : 


f ’ 
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where a z = hF"(0) and b j = H*(0). 

The solutions of the zeroth order equations are now straight- 
forward and yield the following: 

(21a) 

and g o = b x n (21b) 

B. DETERMINATION OF THE HIGHER ORDER COEFFICIENT FUNCTIONS 
The first and higher order coefficient functions f ^ , and 
for the stream function and temperature distribution inner 
expansions are obtained by solving the sets of coupled fifth 



f = a n‘ 
0 2 
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order ordinary differential equations (19) and (20). These 
equations were obtained by equating like powers of E, in the 
momentum and energy expansion defined in equations (17) and 
(18). The presence of e* in the inhomogeneous parts of the 
differential equations introduces another independent para- 
meter which significantly complicates the numerical procedure 
employed in their solution. To eliminate this problem and 
isolate the effects of e* the coefficient functions can be 
redefined in terms of universal functions which depend only 
on the Prandtl number. The revised forms of f^ and gj. appear 
in equations (22), presented below. 
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( 22 ) 



ing only on the Prandtl number. Since the first and higher 
order equations of (19) and (20) ate linear, the principle 
of superposition applies and the equations of (22) can be 
substituted into (17) and (18) and separate like power of 
£* and obtain the following equations. 

The first order momentum equation becomes: 

(23a) 
(23b) 



M [Y ] = - 1 
1 i o 



and M [Y ] = - 4f", 

111 o 

where M is the linear differential operator defined by: 

l 

M = Y'" + 2f Y" - 3f Y' + 3f"Y (23c) 

1 1 0 1 0 1 O' 1 
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The first order energy equation becomes: 



E [Z 3=0 



and E [Z ] = - 4 Z ’ , 

! 11 0 0 

where E is the linear differential operator defined by 
i 

E = i-Z" + 2f Z' - 2 f ' Z + 3Z' Y - Z Y’ 

Pr i oi oi ooi ooi 

The second order momentum equation becomes : 
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where M is the linear differential operator defined by 
2 

M = Y ' ' ’ + 2f Y" - 4 f ' Y ' + 4 f’’Y . 

2 2 02 02 02 

The second order energy equation becomes: 

E [Z ] = 2Y ’ Z - 3Y Z' ; 

2 2 0 1 0 1 0 1 0 1 0 

E [Z ] = 2Y ’ Z + 2Y’ Z - 3Y Z' - 3Y Z' 

2 21 1011 1110 1011 1110 




and E [Z ] = 2Y ' Z - 3Y Z' - 4Z' , 

222 1111 1111 11 

where E is the linear differential operator defined by 
2 

E = i- Z" + 2f Z’ - 3f’Z + 4Z ' Y - Z Y' 

2^ r 2 02 02 22 002 

The following boundary conditions must be satisfied 
equations (23)-(26): 
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Y ' 
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z 



i o 



-*■ 0, Z -> 0; 



(28) 



n 



2 



2 1 



and 



Z 



— -> o , z -> 0, z 

3 2 1 



2 2 - 



which follow directly from the coefficient function boundary 
conditions given in equations (B-33), (B-35) and B-56). 

Although they were not used in the present analysis, the 
third order universal function equations are presented in 
Appendix C. 

C. CONSTRUCTION OF THE OUTER SERIES EXPANSIONS 

The inner series developed in the preceding section is 
useful only in approximating local behavior at the surface in 
the vicinity of the discontinuity. In order to describe the 
temperature and velocity profiles at greater distances from 
the plate another series is required which will satisfy the 
asymptotic boundary conditions and also match the inner series 
at some intermediate distance from the surface. According 
to the method of Goldstein [17] the outer series expansions 
for the stream function and temperature distribution should 
be of the following form: 




(29) 



and 




(30) 
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Since as £-*-0 the original initial conditions must be satisfied, 
the following equalities must hold; thus defining the zeroth 
order coefficients: 



ip (_y*) = F(y*) 

0 



and 



( 31 ) 



0 (y*) = H(y*) 
0 



where F and H are the free convection similarity functions for 
the stream function and temperature distribution, respectively. 
This condition is merely a mathematical restatement of the phy- 
sical condition which requires that the temperature and velocity 
profiles at the discontinuity, £= 0, match those generated by 
the free convection process in the region between the leading 
edge and the discontinuity. 

Now the asymptotic outer series froms ( 29 ) and ( 30 ) are 
substituted into the governing equations (6) and ( 7 ). Once again 
through the process of collecting terms which are coefficients 
of like powers of £, the following set of linear ordinary dif- 
ferential equations is obtained from the momentum equation: 

ip ' ip ' — ’p " ip = 0 ( 32 a) 
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+ 0' - i p 0' - 0' = L_0" (33c) 

2 03 12 2 21 “ 12 21 2 30 rT 2 

The details of the solution of equations (32) and (33) can be 
found in Appendix D and the results are summarized below: 
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With these coefficient functions determined, velocity and 
temperature profiles for large y* can be generated according 
to the following expansions: 



U* - ft. - F'Cy*) * + [§F'(y») 

2 

1 D 

+ ^2 - 2y*) F" (y*) ] £ 3 



(36) 



and 

C 1 D 

0 = H (y *) + ^’(y*)? 2 + - 2y*)H ' (y*) t , 3 (37) 

2 2 

D. NUMERICAL SOLUTION OF UNIVERSAL EQUATIONS 

The universal equations (25) -(26) with boundary conditions 
(27) and (28) were solved numerically with an IBM 3033 series 
digital computer. Double precision machine arithmetic was 
used for all computations. A least squares correction method 
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developed by Nachsteim and Swigert [20] was used as shooting 
point routine to determine the unknown starting values for 
each set of fifth order two point boundary value equations. 

A standard fourth order Runge-Kutta initial value integration 
routine was employed to tabulate values for the free convec- 
tion functions F, F' , F", H and H'; the zeroth order inner 
series coefficient functions f , f', f " , g and g'; and the 
higher order universal functions Y , Y , Y , etc. These 



10 11 2 0 

tabulated values were then used to generate the inner and 
outer velocity and temperature profiles as well as the local 
Nusselt number at the surface of the plate. 
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IV. RESULTS 



A. PRESENTATION OF RESULTS 

1 . Velocity and Temperature Profiles 
a. Inner Series Profiles 

The inner series expansions for velocity and 
temperature were computed form the tabulated values of the 
universal functions which were obtained form the numerical 
solution of the fifth order ordinary differential equations 
defined by equations (23) through (25). The velocity profiles 
were obtained for specified values of e* the blowing or suction 
parameter according to the following equation: 

u* = 5[f’ + (Y' + e *Y ' )? + (Y ' + e *Y ' + e* 2 Y' )£ 2 ] (38) 

0 10 11 20 21 22 

where the f ' , Y' , Y' , etc., are functions of q. In order 
0 10 11 

to facilitate comparison and matching of the inner and outer 
profiles the nondimens ional velocity u* was plotted as a func- 
tion of y* where the relationship between y* and q is: y*=£*q. 
Numerous sets of profiles were generated for analysis and 
comparison; however, only the profiles corresponding to the 
following parameters are presented in this paper: 

1) Blowing and suction parameter, e* , of magni- 
tudes ±0.05 and ±0.25 where positive values cor- 
respond to suction and negative values to blowing. 

2) Streamwise nondimens ional parameter, $, corres- 
ponding to x/L ratios of .25, .50, and 1.0. 
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The plots of u* versus v * for these parameters are presented 
in Figures 2-13 and for ease of comparison both inner and outer 
profiles are included on the same plot. 

The temperature profiles generated by the inner 
series expansions were obtained from the following equation: 

e = l + ?[g + (z + e *z k + (z + c*z + c * 2 z K 2 ] 

0 1 0 1 1 2 0 2 1 2 2 

(39) 

where the g , Z , Z , etc are functions of q; hoi^ever, as 
0 10 11 

was the case with the velocity profiles, 0 was plotted as a 
function of y* with y* = ?*n- The same parameters were used 
in plotting 0 versus y* and the outer profiles were included 
for comparison. These plots appear in Figures 14-25. 
b. Outer Series Profiles 

The outer series expansions which were constructed 
in the preceding chapter approximated the behavior of the ve- 
locity and temperature distributions at greater distances from 
the plate. The outer velocity profiles are given by equation 
(36) which is written below for convenience: 

u* = F ' (y *) + [yL- F" (y *) ] £ 2 + 

[§F’(y*) + \ (^ - 2y*)F"(y*)K 3 (36) 

where F 1 and F" are the first and second derivatives of the 
free convection stream function evaluated at the location of 
the discontinuity, x = 0. They were numerically calculated 
for values of the free convection similarity variable which 
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is identically equal to y* at x = 0. The constant a =»F M (0), 

and the constants C and D are computed numerically from 

1 2 

equations (D-16) and (D-17) respectively. The outer velocity 
profiles corresponding to the same e* and £ parameters as were 
specified for the inner profiles are shown in Figures 2 through 
13. 

The outer expansion which defines the temperature 

profile valid for large values of y* was given in equation (37) 

and is repeated below: 

C , D 

6 = H(y *) + ^i-H'ay*)? 2 + - 2y*)H'(y*K 3 (37) 

2 2 

where the constants C , D , and a , are defined as before, and, 

1 2 2 

the H and H 1 are the free convection temperature distribution 
and its derivative evaluated at x = 0 where v* is the free 

J 

convection similarity variable. Once again the plots of 0 
versus y* for the same parameters as for the inner profiles, 
are presented in Figures 14 through 25. 

2 . Heat Transfer Characteristics at the Surface 

The standard method of characterising the heat transfer 
at the surface of a heated plate is to calculate the local 
Nusselt number at the point of interest. In this particular 
problem the Nusselt number at any location above the disconti- 
nuity, x>0, is given by the following: 

Nu = - _ L — ^ |_ (40) 

(V T J sy y= 0 
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which reduces to the following form when the appropriate sub- 
stitutions are made for T and y in terms of nondimens ional 
variables : 

Nu — Ct 1 ) [g'CO) + (Z' (0) + e *Z ' C0)K 

^ 0 10 11 

+ C Z 1 (0) + £ * Z 1 (0) + £ * 2 Z ’ (0))? 2 ] (41) 

2 0 21 22 

where the universal functions Z* , Z ' , etc. were numerically 

10 11 

evaluated at n = 0 . The heat transfer data is presented in 
graphical form using the following nondimens ional parameter 
Nu 

r cm 

l 4 J 

which was plotted as a function of £ for the following values 
of the suction and blowing parameter: 

£* = ±0.05 and ±0.25. 

The results are presented in Figures 26-29. 

B. DISCUSSION OF RESULTS 

1 . The Matching Condition 

a. The End-Point of Numerical Integration. 

The end-point of numerical integration is the 
value supplied to the Nachsteim and Swigert starting value 
routine at which it attempts to satisfy the boundary conditions 
at n = 00 . The choice of this point has an obvious direct in- 
fluence on the starting values, f"(0) and g'(0), and also on 

k k 

the outer series constants C , C , D , D and D defined 

i 0 1 1 2 0 2 1 2 2 

in ( C - 1 6 ) and (C-17) and which were numerically evaluated at 
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the value of n chosen as the end-point. A summary of the 
effects of varying the end-point on both the starting values 
and outer series constants is presented in Table I. For each 
of the end-point values, H en( j> listed in the table, a complete 
set of tabular profiles were generated and qualitatively com- 
pared to determine the degree of matching evidenced. The 
results of this process revealed that the smaller the value of 
r ’end > * 3etter was t ^ ie overall matching of the profiles. 

However, for n en( j <2.0 no further improvement in the matching 
was observed. Also, the boundary conditions at n = OT were 
satisfied to equal degrees of accuracy for each F enc j selected. 
Therefore, all of the numerical results presented were based 
on an end-point value of F en( j = 2.0. It was also observed 
that the outer profiles were quite insensitive to variations 

in the constants C , C , etc,; and that for a fixed e* and 

10 11 

£ the parameter Nu(Gr/4) was insensitive to changes in 

the initial values which were caused by changing n en< ^- 
b. The Inner-Outer Profile Match Point 

In addition to the qualitative numerical procedure 
outlined above to effect a matching of the inner and outer 
profiles, another approach was undertaken based on the appli- 
cation of an overall energy balance to the problem under con- 
sideration. The principle of conservation of energy requires 
that all energy entering a specified control volume must leave 
that volume providing there is no storage of energy within the 
volume. In the case of the present problem all of the energy 
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TABLE I 



SUMMARY OF STARTING VALUES AND MATCHING CONSTANTS 



'"'end 


2.0 


3.0 


5.0 


7.0 


10.0 


Y" (0) 
1 0 


IO' 9 


10" 9 


IO' 9 


0 


0 


INI 

o 

O 


io " 9 


io " 9 


io ' 9 


0 


0 


Y" (0) 
ir 


2.5787 


2.8611 


3.0893 


3.1745 


3.2263 


Z* (0) 
ip J 


-1.5531 


-1.7623 


-1.9329 


-1.9994 


-2.0403 


Y" (0) 
20 v 


io ' 4 


10' 4 


io ' 4 


io ' 5 


io ' 10 


r (0) 
2 0 v 


io ' 4 


io ' 4 


io ' 5 


io ' 5 


io ' 8 


Y" (0) 
21 V 7 


-2.5437 


-3.3282 


-4.0277 


-4.3080 


-4.4860 


Z' (0) 


.12905 


.21959 


. 345070 


. 4034 


. 44235 


Y” (0) 
2 2 v 


2.15407 


2.3655 


2 . 7645 


2.9743 


3.1229 


Z’ (0) 
2 2 v 


-.96015 


-.9905 


-1.18464 


-1.3222 


-1.4087 


C 

1 0 


0 


0 


0 


0 


0 


C 

l l 


1.605 


1.636 


1.622 


1.615 


1.61 


D 

2 0 


0 


0 


0 


0 


0 


D 

2 1 


3.2534 


4.403 


7.022 


9.875 


14.289 


D 

2 2 


-2.8477 1 


1 -2.9755 


-2.9811 


-2.8945 


-2.788 
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entering the system at the surface of the plate must exit in 
an upward direction parallel to the plate. The mathematical 
formulation of this requirement is presented below in integral 
fo rm. 

00 _ ___ x+L _ _ 

I c P (T - TJudy = f q" (x) dx (42) 

•'0 P *00 

The left hand integral must be divided into two separate in- 
tegrals since there exists an inner and outer region where 
the velocity and temperature behavior are described by differ- 
ent expans ions -- the inner and outer profiles. The dividing 
line between these two regions must, therefore, define the 
match point of the inner and outer profiles. The details of 
evaluating the integral equation (42) are presented in 
Appendix E. The only quantative result of the energy balance 
analysis was that an exact balance could be achieved only if 
the matching point between the inner and outer expansions is 
taken to be y* = 0. That is the integral in (42) must be 
evaluated using the outer series to account for energy leaving 

the volume the left side- -and the inner series to account 

for the thermal energy entering the volume at the surface of 
the plate. As the profiles themselves suggest, the outer 
profiles when viewed in an integrated or global sense more 
nearly satisfy the energy balance, although in a pointwise sense 
they clearly do not, especially at the surface of the plate. 
Since behavior at the surface of the plate whether with respect 
to heat transfer or shear stress is of primary interest, from 
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an engineering standpoint, knowledge of the exact location of 
the match point is not essential, since the inner series ex- 
pansions are always used to furnish results at the surface. 

2. Velocity Profiles 
. - ■ . 

Before analyzing the velocity profiles some method of 
characterizing the relative magnitude of the blowing or suction 
velocity should be established. To accomplish this v^ will be 
compared to the maximum value of the u-component of velocity 
in the case of unperturbed free convection evaluated at the 
location of the discontinuity, x = 0, which will be designated 
by . The ratio v /u which characterises this comparison 

is reduced to the following nondimens ional form through 

appropriate substitutions: 

v ^ 

= 5.124 e*Gr~* (43) 

u 

max 

For a Gr * 10 8 for which the laminar flow regime is still 
applicable the relationship expressed in (43) becomes: 

v 

^ — = .0512 e* (44) 

u 

max 



From this correlation the blowing and suction parameters used 
in the present problem correspond to the following specific 
relationships : 

e* =0.05 implies v = . 002 56 u , 
r w max’ 



and 

e* =0.25 implies v = .0128 u 

r w max 



(45) 
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These values for v are consistent with the original assumption 
underlying the governing boundary layer equations. 

A careful review of the inner profiles plotted in 
Figures 2-13 reveals the following features about velocity be- 
havior with suction or blowing: 

1) As suction is increased at a specified value 
of j-, the location of the maximum u* is drawn 
in closer to the plate and its magnitude is 
decreased. 

2) At a given value of suction the location of 
the maximum u* is moved closer to the plate 
and the magnitude of u* is decreased as the 
value of is increased. 

3) The results for blowing exhibit the reverse 
effects on the location and magnitude of u* 
as those described in 1) and 2) for suction. 

Numerical data on the effects of suction and blowing on 
the location and magnitude of the maximum u* for the inner 
velocity profiles as well the corresponding data for the 
unperturbed free convection case is presented in Table II. 

Comparison of the suction data with the free convection 
data reveals that the location of the maximum u* is shifted 
closer to the plate and also that its magnitude is reduced by 
15% when =0.25 and by 37% when £ = 1.0 with a suction mag- 
nitude of c * =0.05. When the suction is increased to £*=0.25 
the effects are correspondingly amplified. Although this is 
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TABLE II 



LOCATION OF INNER PROFILE MAXIMUM VELOCITIES 





X 

L .25 


. 50 


1.0 


Free Convection: 








y* 


1.020 


1.070 


1.15 


u* 

max 


.3088 


.3383 


. 3907 


Suction c * = 0 . 05 








y* 


. 8356 


. 8221 


. 7814 


u* 

max 


. 2623 


.2568 


.2444 


Suction e *=0 . 25 








y* 


. 6868 


. 6634 


.4906 


u* 

max 


. 2426 


.2194 


.1862 


Blowing e*=-0.05 








y* 


. 8929 


. 9086 


.9267 


u* 

max 


.2665 


.2717 


. 2 851 


Blowing £ *=-0 . 25 








y* 


.9780 


. 9951 


1.036 


u* 

max 


.2589 


.2745 


. 5333 
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the behavior which would be expected for suction, the percent 
reduction in maximum velocity predicted by the inner expansion 
is probably too great considering the relative magnitude of 
suction applied. In the case of blowing the inner series 
fails to produce results consistent with expected physical 
behavior since neither the location of the maximum u* is moved 
farther from the plate nor is the magnitude of u* increased 
above that observed for free convection. Also the degree of 
matching between the inner and outer profiles lessens as the 
suction magnitude or distance from the discontinuity is 
increased . 

As indicated above the results obtained in the case 
of suction are consistent with the physical situation since 
in this case heated fluid is drawn closer to the plate where 
viscous forces tend to reduce the maximum velocity in the 
streamwise direction. In the case of blowing heated fluid 
is pushed farther from the surface away from the influence 
of viscous forces and into the region where buoyancy forces 
act to accelerate the flow in the streamwise direction there- 
by increasing the maximum u*. In addition for a specific 
value of either suction or blowing the effects described 
above will be accentuated as the distance from the disconti- 
nuity is increased; i.e.; for increasing values of j- or 
since the blowing or suction effects will have had a greater 
opportunity to predominate over unperturbed free convection. 
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The lack of positive correlation in the case of blowing versus 
unperturbed free convection might be attributed to any of the 
following : 

1) Insufficient terms in the inner series expan- 
sion; since the signs of the coefficient 
functions alternate in a pairwise manner, 
addition of another term would possible serve 
to effect a better correlation. 

2) The effects of blowing become significant in 
the region where the inner series no longer 
provides a valid representation of the velocity 
behavior . 

The effects of suction and blowing on the outer velocity 
profiles are graphically presented in Figures 2-13; however, 
the specific effects on the location and magnitude of the 
maximum value of u* are given in Table III using the same 
format as Table II for the inner profiles. The results in 
the case of suction are somewhat inconsistent with the behavior 
expected relative to the free convection case. The location 
of the maximum velocity does not move closer to the plate, 
and its magnitude does not decrease when compared to the - free 
convection profiles. This anomaly suggests that the effects 
of suction on velocity are more significant in the region 
close to the surface of the plate where the behavior is better 
approximated by the inner series expansion. Conversely, in 
the blowing case the results correlate well with expected 
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TABLE III 



LOCATION OF OUTER PROFILE MAXIMUM VELOCITIES 





X 

L’ . 25 


. 50 


1.0 


Free Convection: 








y* 


1.020 


1.070 


1.150 


u* 

max 


. 3088 


. 3383 


. 3907 


Suction £*=0.05 








V* 


1.053 


1.1588 


1.3628 


u* 

max 


. 3110 


. 3466 


.4212 


Suction e*=0.25 








y* 


.9586 


. 9519 


1.1072 


u* 

max 


. 3117 


. 3468 


. 4160 


Blowing £*=-0.05 

y* 


1.110 


1.240 


1 .472 


u* 

max 


. 312 


. 350 


.433 


Blowing £*=-0.25 

y* 


1.206 


1.370 


1.617 


u* 

max 


. 317 


. 364 


. 466 
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Figure 2. Suction Velocity Profiles: e* = 0.05, = 0.25 
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Figure 3. Suction Velocity Profiles: e* = 0.05, ^ = 0.50 
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Figure 4: Suction Velocity Profiles: e* = 0.05, 
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Figure 5: Suction Velocity Profiles: e* = 0.25, 
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Figure 6: Suction Velocity Profiles: c* = 0.25, 
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Figure 7: Suction Velocity Profiles: e* = 0.25, 
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Figure 8: Blowing Velocity Profiles: e* = -0.05, 
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Figure 9: Blowing Velocity Profiles: e* = -0.05, 
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Figure 10: Blowing Velocity Profiles: e* = -0.05 
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Figure 11: Blowing Velocity Profiles: e* = -0.25, 
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Figure 12: Blowing Velocity Profiles: e* = -0.25, 
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Figure 13: Blowing Velocity Profiles: c* = -0.25, 
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behavior. Relative to the free convection case the maximum 
value of u* moves away from the plate, and its magnitude in- 
creases from 1.01 to 10. 8%, with a blowing parameter e*=-.05, 



as — increases from .25 to 1.0. When the blowing magnitude 
Li 

is e* = -0.25, the value of u* 



max 



increases from 2.7% to 19.2% 



x 



as increases from .25 to 1.0. Whether or not this increase 
is excessive must await comparison with experimentally obtained 
data. In any case it appears that a very small perturbation 
of the momentum boundary layer at the wall, produces relatively 
significant changes at increased distances from the plate. 

3 . Temperature Profiles 

The effects of blowing and suction on the temperature 
profiles are not as visible in the plots Figures 14-25, as 
was the case for velocity since there are no maximum velocity 
peaks to observe. However, the same general pattern emerges 
as did with the velocity profiles, and key points are summarized 
below for the inner temperature profiles: 

1) For a given location above the discontinuity the 
inner profile pivots clockwise about 8 = 1.0 as 
the magnitude of suction is increased, this re- 
presents an increase in the gradient which is 
characteristic of an increased heat transfer rate. 

2) For given magnitude of suction there is a marked 
increase in gradient as the distance from the 
discontinuity in the streamwise direction is in- 
creased, i.e., the value of ^ increases. 

Li 



61 




62 



ti x| 



t 




Figure 15: Suction Temperature Profiles: e* = 0.05, 
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Figure 16: Suction Temperature Profiles.: e* = Q.0.5, 
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Figure 17: Suction Temperature Profiles: e* = 0.25, 
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Figure 18: Suction Temperature Profiles: e* = 0.25, 
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Figure 19: Suction Temperature Profiles: e* = 0.25, 
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Figure 20: Blowing Temperature Profiles: e* - -0.05, 
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Figure 21: Blowing Temperature Profiles: e* = -0.05, 
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Figure 22: Blowing Temperature Profiles: e* = -0.05, 
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Figure 23: Blowing Temperature Profiles: e* = -0.25, 
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Figure 24: Blowing Temperature Profiles: e* = -0.25, 
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Figure 25: Blowing Temperature Profiles: e* = -0.25, 
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Relative to the unperturbed free convection profiles 

•>£ 

the effect of suction for a given value of is to 
increase the gradient. 

4) For blowing the reverse effects of those outlined 
in 1) through 3) above are observed with the addi- 
tional feature that the inner profiles are bowed 
out as well decreasing their gradient. This be- 
havior characterizes a decreased heat transfer 
rate . 

These observed phenomenon correspond with what physical 
intuition suggests should occur. In the case of suction heated 
fluid is moved closer to the wall thereby reducing the thick- 
ness of the thermal boundary - layer and enhancing the rate of 
heat transfer as is evidenced by the increase in gradient. The 
opposite situation in which blowing forces heated fluid 
farther from the plate causes the thermal boundary - layer to 
thicken. This creates an insulating effect which serves to 
reduce the rate of heat transfer evidenced by a decrease in 
gradient . 

The outer temperature profiles again reflect the 
distinctive shape of unperturbed free convection profiles. 

In the case of suction their behavior is inconsistent and 
they do not behave relative to the free convection profiles 
as do the inner profiles. This again suggests that the effects 
of suction are most evident near the surface of the plate 
where the inner series is a closer approximation of the actual 
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behavior. For the case of blowing the outer profiles consis- 
tently move away from the plate with increases in either 
magnitude or distance above the discontinuity. Also, rela- 
tive to the free convection profiles, they behave as expected 
moving farther out from the surface as the magnitude of 
blowing is increased for a specified value of This suggests 

that the outer profiles are indeed a good approximation of 
temperature behavior at large distances from the plate. 

4 . Heat Transfer at the Surface of the Plate 

In order to analyze the effects of discontinuous suc- 
tion and blowing on the rate of heat transfer at the plate a 
graphical comparison was made of the following three physical 
situations : 

1) Unperturbed free convection from an isothermal 
plate . 

2) Free convection from an isothermal plate with 
continuous suction or blowing applied over the 
entire length, (from the study by Sparrow and 
Cess [ 7] . 

3) Free convection from an isothermal plate with 
discontinuous suction or blowing applied at a 
distance L above the leading edge. 

For each of these cases the heat transfer data was correlated 
according to following nondimens ional parameter: 

Nu (Gr/4) ' % (46) 
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where the Nusselt number in each case Nu , Nu, and Nu was 

o i 

evaluated according to length scale associated with the 
discontinuity; i.e., L. The three parameters have the follow- 
ing form: 

1) Nu (Gr/4) = - H ' (0) (x + l)'* 5 (47) 

0 

2> Nu (Gr/4)'^ = - H’(0)(x + 1) - 1.40e* (48) 

1 

3) Nu(Gr/4)'^ = - [H’(0) + (Z’ (0) + e*Z’ (0)5 

10 11 

+ (Z’ (0) + e *Z ' (0) + e * 2 Z ' (0) ) C 2 ] (49) 

20 21 22 

The heat transfer parameters were plotted as functions of x, 
i.e., j-, for suction and blowing magnitudes, e*, of .05 and 
.25 in Figures (26) -(29). As expected in the region from the 
leading edge, x = - L, to the discontinuity at x = 0 the heat 
transfer rates for cases 1) and 3) are identical and in case 
2) the rate is higher by a fixed amount over the unperturbed 
free convection case. At the discontinuity the curve for the 
discontinuous case experiences an abrupt change in slope and 
intersects the continuous case curve at approximately = 

0.22 for e * = 0.05 and at ^ = 0.52 for £* = 0.25. The magni- 
tude of applied suction apparently has a significant effect 
on the location of the crossover point where the discontinuous 
heat transfer rate exceeds that for the continuous case. The 
behavior of the discontinuous heat transfer rate as t- increases 
does not appear to approach an asymptotic suction limit charac- 
teristic of a constant boundary-layer thickness. However, 
since the inner temperature series expansion constructed to 
approximate behavior near the surface is only valid in the 
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Figure 26: Comparison of Heat Transfer at the Surface of the 

Plate with Suction, e* = 0.05. 
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Figure 2 7: Comparison of Heat Transfer at the Surface of the 

Plate with Suction, e* = 0.2S. 
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Figure 28: Comparison of Heat Transfer at the Surface of the 

Plate with Blowing, e* = -0.05. 



79 



1.0 




Figure 29: Comparison of Heat Transfer at the Surface of the 

Plate with Blowing, e* =-0.25. 
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neighborhood of the discontinuity, i.e,, small values of 
it cannot be expected to accurately predict results in asymp- 
totic limiting case as j- becomes large. The addition of higher 
terms to the inner expansion would serve to extend the range 
of validityj however, since the Sparrow and Cess analysis of 
the continuous case employed only first order terms, it is 
reasonable to expect the results of the present analysis to 
provide a closer approximation to the physical situation. 

This essential difference between the two studies should be 
kept in mind when comparing results. The situation with 
blowing applied as depicted in Figures (28) and (29) shows the 
expected reduction in heat transfer rates when compared to 
the unperturbed free convection in both the continuous and 
discontinuous cases. At the discontinuity, however, the curve 
representing discontinuous blowing exhibits a downward shift 
in slope and then follows a path which is basically parallel 
to the free convection and continuous blowing curves. The 
crossover that occurs in the case of a blowing magnitude of 
e* = -0.05, is probably an anomaly due to an insufficient 
number of terms in both the continuous and discontinuous series 
expansions. For a blowing parameter of c* = -0.25 all three 
curves appear to asymptotically approach the same slope. A 
final judgement as to the degree of accuracy with which the 
present analysis models the heat transfer behavior with discon- 
tinuous blowing or suction applied must await comparison with 
experimental data. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



An overall review of the results obtained from this analysis 
leads to a number of general conclusions. An apparent range of 
validity for the expansions can be estimated from a comparison 
of the inner and outer profiles presented in Figures 2-25 based 
on the degree of correlation between them. In terms of the 
streamwise distance from the discontinuity a qualitative upper 

X 

bound of j- = 1.0 is suggested. In terms of the blowing and 
suction parameter the relationship given on equation (44) sug- 
gests an upper limit of |e*| = 1.0 which equates to a magnitude 
of approximately .05 u for v for which the original boundary 
layer equations should remain valid. In the case of suction 
the inner profiles appeared to model the expected behavior of 
the location and magnitude of the maximum velocity with respect 
to the unperturbed free convection profiles better than did 
the outer profiles. When blowing was applied, the outer pro- 
files correlated with anticipated behavior in this respcet 
better than the inner profiles. In general the matching of 
the inner and outer velocity profiles improves as the stream- 
wise distance from the discontinuity decreases and is also 
qualitatively better in the blowing cases. The analysis also 
indicated that while the inner series provides reasonably good 
results at very small distances from the plate, in an overall 
integrated sense the outer profiles exhibit a closer correlation 
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with the velocity behavior associated with free convection. 

The effects of suction and blowing on the temperature 
distribution both appear to be well modeled by the inner 
series expansions. Again the degree of matching between the 
inner and outer profiles significantly improves as the stream- 
wise distance from the discontinuity decreases. Also the 
effects of the discontinuity on temperature behavior, relative 
to the free convection case, are less direct than on velocity 
behavior. This difference is evidenced in the manner of coup- 
ling exhibited by the ordinary differential equations which 
define the inner expansion coefficient functions (see equations 
(19) and (20)). This fact in addition to the apparent over- 
estimate of velocity magnitudes produced by the outer series; 
the underestimate by the inner series; and the heat transfer 
behavior predicted for small suction magnitudes suggest that 
additional higher order terms be added to both inner and outer 
expansions to achieve a better approximation of the physical 
situation. In general the results obtained for the effects 
of discontinuous suction or blowing on heat transfer at the 
surface of the plate seemed consistent with the results for 
unperturbed free convection and the case of continuous suction 
or blowing applied over the entire plate length. Comparison 
with experimental data when it becomes available will enable 
a more definitive conclusion to be made regarding the validity 
of this solution with respect to heat transfer performance. 



83 



Also, with the addition of higher order terms to the 
inner series it would be possible to determine if the heat 
transfer results predicted by the present analysis tend to 
approach the asymptotic limits for large values of £ es- 
tablished by Merkin [9] . In that study the following 
asymptotic limits for continuous suction and blowing at 
large distance from the leading edge were reported, given 
below in terms of the parameters used in this analysis: 

1) For suction Nu/(^— -»■ 4e*Pr as ^ -> 00 

2) For blowing Nu/ (j^-) 4 ->■ 0 as -*■ 00 

The asymptotic limit in the case of suction is represented 
by the dashed line in Figures (26) and (27) . 

Although an exact determination of the inner-outer series 
match point is not essential when the primary interest of 
the investigation is to predict behavior at the surface of 
the plate, a more rigorous determination of that point would 
serve to enhance the overall creditability of the present 
method of solution. Therefore, another recommendation to 
improve the method point is to use the energy balance 
integral given in (42) to form an energy difference, i.e., 
from a difference integral by substracting the right side 
from the left side. A numerical minimization procedure 
could then be applied to this difference to determine the 
match point which most closely satisfies the energy balance. 
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APPENDIX A 



CONSTRUCTION OF THE FREE CONVECTION MATCHING SERIES 
The equations which define the behavior of the free con- 
vection velocity and temperature profiles were presented in 
equations (9a) and (9b), with appropriate boundary conditions, 
in terms of similarity functions F and H which characterize 
the stream function and temperature distribution respectively. 
In the region near the surface of the plate, y* = 0 , and at 
the location of the discontinuity, x = 0, these functions can 
be approximated by suitable series expansions in powers of 
y* as follows : 



where the variable y* can be used since it is identically 
equal to the free convection similarity variable at x = 0 . 

The constant coefficients in (A-l) and (A-2) can be 
determined from applying the following previously enumerated 
boundary conditions: 

At y* = 0: F(0) = 0 



F(y*) = S a.fy*) 3 
j=o J 



(A-l) 



H (y*) =1+2 b.(y*)j 




(A-2) 



F' (0) = 0 



H (0) = 1 



(A-3) 



As y* ■+ 00 : F' (°°) = 0 



H (°°) = 0 
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These conditions immediately yield the result that a = a = 0. 

o 1 

By evaluating the series representations for F" and H' the 

constants a and b are specified as: 

2 1 

a = (0) and b = H' (0) 

2 1 

where both of these initial values are well known and depend 
only the Prandtl number. 

To evaluate the higher order coefficients the form for 
F''' must be obtained from equation (9a): 

F' ' ' = - 3FF" + 2F ' 2 - H 

which when evaluated at y* = 0 yields F'''(0) = - 1 and 
equating this value to the series form of F 1 '' yields: 



A similar technique is used in determining b from equation 

2 

(9b) and the following value is obtained: 



The remainder of the coefficients are found by successively 
differentiating equations (9a) and (9b) and equating them to 
the appropriate series representation evaluated at y* = 0. 
The results of this process are summarized below: 



1 




b = 0 



2 



a 



0 



b = H' (0) 



0 



a 



0 



b = 0 



2 



(A- 4) 



a = lsF M (0) 
2 



b 



3 



0 



a 



3 



1/6 



b 



4 



-PrF" (0) H ' (0) 
8 
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a 

4 


-H* CO) 
24 


a 

5 


[ F" ( 0 ) ] 2 
120 



These series forms for F and H will be used to satisfy the 
matching conditions for both the inner and outer expansions 
later in the course of the analysis. 
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APPENDIX B 



DETERMINATION OF THE STREAM FUNCTION AND TEMPERATURE DISTRI- 
BUTION SERIES EXPANSION PARAMETERS 

Recalling that the definitions of the stream function and 
temperature profile series expansions involved certain arbi- 
trary constants, the methods and procedures used in evaluating 
them are outlined in this appendix. 



1 



= e*(?/<J> ) a + * C c U k f f (n)] 

i 3 r 


(B-l) 


9 = 1 + <D 4 C d [? k g k (n)] 


(B - 2 ) 


where 




o K 

£ = d> x and q = d> y *x 
1 2 


(B-3) 



The object of constructing such expansions for the stream 
function and temperature distribution is to obtain series re- 
presentations which satisfy the momentum and energy equations 
which mathematically model the physical problem. The indivi- 
dual components of those governing equations (6) and (7), i.e., 
the partial derivatives, must be evaluated in terms of 
expressions (B-l) through (B-3). Substituting ( B - 1 ) and (B-2) 
into the momentum and energy equations yields: 

Momentum Equation 

1 -2b 2b - 1 

4$ a + a { [a(c + k) + b)£ k f£] [£ k f£] 
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_2b 

[a (c + k)? k f k ] [C k f^] } - 4e*0 j a ^ 



c + 



a u k fj;] 



2b 



_3b c + — 

4 , a <j) 3 b C a [? k £, " ' ] + 1 + 4> [C k g k ] 

12 3 K 4 K 



(B - 4 ) 



Energy Equation 




- [a (c + k)5 k f,] [C k g' ]} - 4e *<}) a <{, <j> C 

K K 12 4 



b 




[5 k gil 



k 




2b 2b 



(B - 5) 



Certain conditions which must be satisfied will impose 
restrictions on the allowable values for the <j> ' s and constants 
a, b, c and d. The one condition which must be satisfied by 
the equations is the initial velocity and temperature matching 
at a distance L from the leading edge of the plate where the 
discontinuity in normal velocity begins. At this location the 
perturbed velocity and temperature profiles must match the 
free convection profiles. In terms of the stream function the 
velocity is given by: 



The free convection velocity profile written in terms of u* 
can be expressed as follows: 




ca ac 1 ka ka 



= F' ly *] = Z ja (y * ) •* " 1 

j-2 J 



(B-7) 
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Equating the leading term in each series 



when j=2 and k=0 we obtain the folloxving 

ca ac , 

o „* ,c, b . r y*. b f' 

2a y* = H 4> (t-) 0 

2 12 3b 

For proper matching, the exponents of y* 



ac 

'b“ 



-1 = 1 



(B- 6) and (B- 7) i . e . , 
relationship : 

(B - 8) 

must be equal; therefore, 
(B-9) 



Also, since it is desirable to have expressions which are in- 
dependent of e* for the boundary conditions on the f^ as ir >co , 

the following condition will be imposed; 
ca 



<}) C <f> b <J> = 1 

1 2 3 



(B-10) 



The requirements for matching the temperature profiles yields; 

k 



ad ka ka 

b , r y*. r xk , "b /-y^ b 



b Cj-) b s k ] - jS 1 b j y‘ J (B-ii) 
Applying similar arguments to those used for the condition of 
velocity matching, the following relationships are obtained 



from (B- 11) : 

ad 
"b 



= 1 



(B-12) 



and 



4> d 4> <f> =1 C B - 1 3) 

1 2 4 

For additional restrictions on the exponent constants consider 
equation (B-4) in which the highest order derivative term must 
appear in the Oth order expansion in order to obtaine meaningful 
results. This requirement imposes the following condition on 
the exponents of 



c + 



3b 

a 



< 2c + 



2b - 1 
a 
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or equivalently, 
ac > b + 1 



(B- 14) 



The same restriction also applies to the energy expansion, 
equation (B-5) and the same relationship is obtained (B-14). 
Since fractional exponents on E, are not permissable, due to 
the inability to match coefficients of like of powers of E, in 
both expansions, from equation (B-12) the following further 
information is obtained: 



Now combining (B-9) and (B-16), the constant c is specified 
as : 



Using this information and the relationship given in (B-14) 
imposes the following restrictions on the possible values of 
a and b : 



The relationships previously obtained for the ' s can be sim- 
plified by substituting the known values for the exponents in 
equations (B-10) and (B-13) from which are obtained the 
following : 



d = 1 



(B- 15) 



and 




(B-16) 



c = 2 



CB-17) 



a > 1/3 and b < - 1/3 



(B - 18) 



4> 2 <j> 2 cJ> = 1 



(B- 19) 



1 2 3 



and 



4> • <j> (p ■ = 1 



(B- 20) 



1 2 ^ 
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Since equation (B-5) is homogeneous in <p and <J> appears only 
singly in equation (B-4) , without loss of generality it can 
be set equal to unity. Setting <j> = 1 in equations (B-19) and 

4 

(b -20) yields the following further relationships for the <p's: 
=1 (B- 21) 

1 2 

and 

<J> = 1 (B- 22) 

3 

The foregoing analysis has permitted the reformulation of 
equations (B-4) and (B-5) in terms of only two unspecified 
parameters a and <t> , and those equations reduce to the follow- 
ing forms: , , 

— 3 - — 

Momentum: 4a<p^ a £ a {[k + 1) £ k f£] [ C k f^l - [ (k + 2)£ k f k ] . 

[C k f^] > - 4c*£[5 k f£] = [C k f^'l + 5 + 5 2 [C k g k ] (B - 2 3) 

1 3 _1 

Energy: 43^ a £ a {[(k + 1) C k g k l U k f£] - [ (k + 2)^ k f k ] * 

[C k g^]) - 4c*£[S k g£] = ^rtC k g”] (B- 24) 

Equation (B-18) allows for only two reasonable choices for 
the value of a; namely, a = 1/3 or a = 1/2. Furthermore, for 
each of these two values <p can be defined to either include 

i 

e * or not; thereby posing four possible cases for consideration. 
The cases which specify a = 1/2 do not satisfy the boundary 
condition for f as and therefore, are unacceptable for 

i 

the present problem. The case for which e* is included in <j> 
and a = 1/3 results in a set of ordinary differential equations 
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which include c* in the homogeneous part and which cannot be 
transformed into general equations independent of z* through 
the use of universal functions. Therefore, the case with 

a = 1/3 and c* not included in the definition of p was 

1 

selected to complete the solution of the problem. A summary 
of the results is presented below: 



a = 


1/3, 


b = 


- 1/ 


3 , c= 2 


, d = 1 






(B- 


25) 






1/3 




S 


1/3 










♦i 


■ # 


9 


<f> = 

2 


CfJ 


, <P 

3 


= <P = 

4 


1 


(B- 


26) 


c = 


(£*) 


1/3 

9 


n = 


(|x) 


1/3 

y* = 






(B- 


27) 


p = 


- e* 
3 £ 


c 3 + 


C 2 K 


kf ki 








(B- 


28) 


0 = 


1 + 


C[C k g k ] 










(B- 


29) 


The momentum 


and 


energy exp 


ans ions 


which 


result from 


these 




def init 


ions 


were 


pres 


ented 


in Chapter I I , 


equations 


(17) 





and (18) . All that remains is to specify the appropriate 
boundary conditions for the f^. and g^.. 

At n = 0: u* = ||* = SU k f£] = 0, (B - 30) 

V * = "H = ‘ [e * + J C" 1 [ Ck+2) ^ k £ k ] } = ±e*, (B- 31) 

and 

9 = 1 + C[? k g k ] = 1 CB- 32) 

From these original boundary conditions, the local boundary 
conditions for f^. and f£ are determined to be: 

f k (0) = f^ CO) = g k (0) = 0, for all k (B-33) 

The boundary conditions at n -> ' CD are obtained through matching 
the velocity and temperature profiles to the free convection 



93 



profiles. The velocity match is performed by equating equations 
(B-6) and (B-7) with the appropriate values substituted for the 
unspecified parameters. The result of this process is the 
following relationship: 



js 2 jy?*) 1 " 1 ■ [O k+lf i' 



CB- 34) 



By equating like powers of y* and substituting the values for 
the a^ obtained in Appendix A the following local boundary con- 
ditions for the f£ as tt * 00 are determined: 



—2. = 2a = F" (0 ) as r\- 

n 1 



f ' 

n 



A = 3a = - i 

2 3^ 



as m°° 



(B - 35) 



— = 4a = - 

n 3 



H' (0) 



as n-* 00 



-2- = 5a = 



[F"(0)] 



24 



as p-*- 00 



The temperature profile matching yields the following local 
boundary conditions for the g^.: 

g„ 

— 2- = b = H ' (0) as tv* 00 
n i 



g, 

= b = 0 
n z 2 



4 = b = o 
n 3 3 

ii. = b = - 

n 4 4 



as 



as n-*-°° 



(B- 36) 



PrF" (0)H ' (0) 
8 



as rv 4,00 
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APPENDIX C 



THIRD ORDER UNIVERSAL FUNCTION EQUATIONS 
The third order ordinary differential equations which 
result from the expansions given in (17) and (18) are pre- 
sented below, although only functions of second order 
and lower ivere employed in the actual analysis: 

f ' ' ' + 2f f" - 5 f ' f' + 5 f" f = 5 f ' f ' - 3f'f" - 4f"f 

3 03 03 03 '12 12 12 

- 4e*f" - g ( C 

2 1 



— g" + 2f g 1 - 4 f ' g + 5g'f - g f' = 2g f' + 3g f' 

P r 3 03 03 03 03 13 21 



- 3f g'- 4f g' - 4e *g ' 
12 2 1 



(C-2) 



The corresponding forms for f and g in terms of the univer- 

3 3 

sal functions and the resulting differential equations are as 
follows : 

(C- 3a) 



f = Y + e *Y + c * 2 Y + £ * 3 Y 

3 30 31 32 33 



and 



= Z + £ *Z + £* 2 Z + £* 3 Z 

3 30 31 32 33 



(C- 3b) 



The third order momentum equation becomes: 

M [Y ] = 5Y ' Y' - 3Y Y" - 4Y" Y - Z ; (C-4a) 

3 30 1020 1020 1020 10 

M [Y ] = 5Y ' Y ' + 5Y ' Y ' - 3Y Y" - 3Y Y" 

3 31 1120 1021 1021 1120 
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(C- 4b) 



- 4Y" Y - 4Y" Y" - Z ; 
1 0 2 1 1 1 2 0 1 1 



M [Y ] = 5Y ' Y ' + 5Y ' Y' - 3Y Y" - 3Y Y" 

3 32 1022 1121 1022 1121 



- 4Y" Y - 4Y M Y - 4Y" 

1 0 2 2 1 1 2 1 2 1 



(C-4c) 



M [Y ] = 5Y ' Y ' - 3Y Y" - 4Y" Y - 4Y" (C-4d) 

3 33 1122 1122 1122 22 

where the linear differential operator M is defined by: 

3 



M = Y"' + 2f Y" - 5f ' v ' + 5 f "y 

3 3 03 O' 3 O' 3 



(C-4e) 



The third order energy equation becomes : 



E 


[Z 


] = 2Z 


Y 1 + 


3Z 


Y' 


- 


3Y 


Z' - 4Y 


Z ' 




3 




3 0 1 


0 2 0 


2 


0 1 


0 


1 0 


2 0 2 


0 1 


0 


E 


[Z 


] = 11 


Y" + 


11 


Y' 


+ 


3Z 


Y" + 51 


Y 


t 


3 




3 1 1 


0 2 1 


1 


1 2 


0 


2 


0 1 1 


2 1 


1 0 






- 3Y Z ' 


- 3Y 


Z' 


_ 


4Y 


1 ' 


- 4Y Z ' 


_ 


4Z T 






10 2 1 


i 


1 2 0 




2 


0 1 1 


2 1 1 


0 


2 0 


E 


[Z 


] = 11 


Y T + 


11 


Y' 


+ 


51 


Y' + 3Z 


Y 


t 


3 




3 2 1 


0 2 2 


1 


1 2 


1 


2 


l i i 


2 2 


1 0 






- 3Y Z ' 


- 3Y 


Z' 


. 


4Y 


Z' 


- 4Y Z ' 


_ 


4Z 1 






10 2 2 


i 


1 2 1 




2 


1 1 1 


2 2 1 


0 


2 1 



(C- 5a) 



E [Z ] = 2Z Y ’ + 3Z Y ' - 3Y Z' - 4Y Z' - 4Z * , 

3 33 1122 2211 1122 2211 22 



(C-5d) 



where the linear differential operator E is defined by: 

3 

E = i-Z” + 2f Z' - 4f'Z + 5 Z * Y - Z Y' . (C-5e) 

3 r r 3 03 03 003 003 

The boundary conditions for these equations are: 



at n = 0 : Y = Y = Y = Y =0; 

30 31 32 33 



and 



Y ' = Y' = Y* = Y' = 0 ; 

3 0 3 1 3 2 3 3 



Z = Z = Z Z =0; 

3 0 3 1 3 2 3 3 



C C - 6 ) 
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and as n^ 00 : 



Y' 



3 Q 



[F"(Q)3 

24 



; Y* + 0 ; Y 1 0; Y' 

31 32 33 



(C-7) 



and lil -> -PrF'COJ.HJJOj z . o; Z - 0; Z + 0 

n 4 8 31 32 33 

which follow from the coefficient function boundary conditions 



given in equations (B-33), (B-35) and ( B - 36) 
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APPENDIX D 



DETERMINATION OF THE OUTER SERIES COEFFICIENTS 
The process of determining the coefficient functions of 
the outer series expansions defined in equations (29) and 
(30) involves the solution of the ordinary linear differential 
equations (32) ; the evaluation of the resulting constants of 
integration; and the solution of the algebraic equations (33). 
The procedure is greatly simplified if the outer expansions 
for ip and 9 are first transformed into series with coefficients 
which are powers series in y* , the coefficients of which are 
derived from the asymptotic forms of the f^. and g^.. This tech- 
nique will facilitate the evaluation of the constants of in- 
tegration significantly reduce the complexity of equations 
(39). 

The asymptotic forms of the inner series coefficients can 
be written as polynomial in n with as yet undetermined 
coefficients as follows: 



f 

2 


= a n 2 
2 




(D- la) 


f 

1 


= A n 3 
1 


+ Bn 2 + Cn + D 

i i i 


(D- lb) 


f 

2 


= A n 4 
2 


+ B n 3 + c n 2+ D n + E 

2 2 2 
o 


(D-lc) 


g 

0 


= b n 

i 




(D- Id) 


g ! 


= B*n + 

i 


c* 

i 


(D-le) 


g 2 


= B*n + 

2 


C*n + D 
2 2 


(D-lf) 
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which must satisfy the asymptotic boundary conditions specified 
in (B-35) and (B-36). The constants are evaluated by substi- 
tuting the equations (D-l) into the original differential 



equations (19) and (20) and 
tions which appear as coeff 
results of this process are 



A 


-i 


A 


. -bi 


i 


6 


2 


24 


B 

i 


= 0 


B 

i 


= 0 


C 

i 


= TBD 


C 

2 


= -Ci 
4a 

2 


D 

i 


= "- e* 
3 e 


D 

2 


= TBD 




then solving the algebraic equa- 

cients of like power of n* The 

summarized below: 

B* = 0 
i 

c * = (D- 2) 

1 2a v ' 

2 

B* = 0 

2 

nit — ~ u i C i 

2 8a 2 

2 



2 6a 2 



where a = i F”(0), b = H'(0) and C and D to be determined 
2 £ 1 12 

from application of the matching condition. This process 
will be outlined at the end of this appendix. 

Now the equations (D-l) with the coefficients as defined 
in (D-2) can be substituted into stream function (15a) and 
temperature distribution (15b) . Then if the definition of 
n = y*/£ is invoked and all coefficients of like powers of £ 
grouped together, the following expansions are obtained ivhich 
are valid as v* becomes large: 

<l> = [a 2 y* 2 - 5V* ! - y* 4 ] * [OH 

* [C y* • & y* 2 H 2 * [D y*]£ J CD-3) 

1 H-d-2 2 
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and 0 = [1 + by*] + [ 0 ] C + 



b D 
1 2 



b C 

i i 
2a 



b C 

_2 L y * 

8a 2 



6a 

2 -J 



(D-4) 



With these equations in hand the final determination of the 
outer series coefficients can now proceed. 

The first of the differential equations defined in (32a) 
is logarithmic and its solution has the form: 

ip = K >p ' CD- 5) 

1 i o 

Since ip is the coefficient of the linear term in the outer 
i 

expansion and that coefficient in (D-3) is zero, K must also 
equal zero, and has the final form: 



ip = 0 
i 



(D- 6) 

Using this equation in (D-2b) reduces it to logarithmic form 



with solution for ip as before: 

2 

ip = K ip ' CD- 7) 

2 2 0 

where K is evaluated by equating (D-7) with the appropriate 
2 

form of the £ 2 coefficient in (D-3). The final result for 

ip then becomes: 

2 C 

•p = — 1— ip ' 

2 2a 0 



CD- 8) 



The last of the differential equations (32c) can be expressed 
in a simpler form using (D-6) ; however, its solution requires 
a different technique. The reduced form for (32c) is: 

(D-9) 



i ip ’ip ' - \ ip " = 0 + ip " ' 

4 0 3 L 



3 0 



100 



Using the free convection stream function equation (9a) (j; " ' 

o 

can be written as: 



if>" ' = - 3ip ip" + Zip ' 2 - 0 
0 0 0 0 0 

and by suitable algebraic manipulations in 

and (D-10) the following differential form 



= 6 



o_ 

ip ' 



- 2 



(D-10) 

combining (D-9) 
is obtained: 

(D- 11) 



\ 0 / ' 0 ' 

Straight forward integration yields the following general 
solution for ip : 



3 



^=6^-2y*^ , +K^' (D-12) 

3 0 0 3 0 

where K is obtained as before with the aid of equation (D-3). 

Thus the particular solution ip has the following form: 

D 3 

ip = 6ip - 2y*ip ' + ^ 5 — ip ’ (D-13) 

3 o • o za o 

2 

This completes the determination of the outer expansions 

for ip and 0 with the exception of the constants C and D . 

1 2 

Since these constants appeared originally in the asmvptotic 

forms of f and f and are determined from the matching of the 
1 2 

inner and outer series, the results of equating the universal 
function forms to the asmyptotic forms will yield expressions 



for C and D 


. The 


following 


equations result 


from this 


1 2 






process : 










For f : 


1 3 

-7- n 


4 

+ C n - T s 


* = Y + e*Y 


(D- 14) 


i 


6 o 


1 0 O 


10 11 


and for f : 


-Ll n 4 


- ^ v ♦ 


C 2 

d n + = 




2 


24 0 


4a 0 


2 0 4a 








2 


2 








Y + e*Y + £* 2 Y 

20 21 22 


(D-1S) 



101 



(D- 16) 



Solving for C from (D-14) yields the following: 

C = C + e*C 
i 10 11 



where 



Y 



C = HH 1 + z- , 



1 o 



and 



Y 



C = (-^ + i-o 
ii n 3n 



A similar procedure yields the following form for D obtained 

2 

from ( D - 1 5 ) : 



D = D +e*D +e* 2 D 
2 20 21 22 



(D-17) 



where 

D 

D 



2 0 



2 1 



Y 

= ( 2 0 
L n 

Y 

= (-— 1 - 



+ - 



1 - 3 , 10 n . 

U n %~a~ n 



1 o 



4a n 
2 



) , 



n 



1 1 
4a 



n - 



c c 

13 1 1 > . 

2a n J ’ 
2 



and 



D 

2 2 



Y C 

r 2 2 . 1 1 

1 n 4a n 



2 

where the value n was taken as n en( j the chosen end point of 

numerical integration. This completes the derivation of the 

outer series expansions, except to note that the constants 

C , C , D , D , and D are evluated numerically from the 
1 0 1 1 2 0 2 1 2 2 

tabulated universal functions. 
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APPENDIX E 



INTEGRAL ENERGY BALANCE ANALYSIS 
In only heat transfer problem an overall thermodynamic 
energy balance must be satisfied. For the case of the present 
problem the integral energy balance takes the following form: 



A p(T - TJ udy = / X+L g" (x)dx 
o P 00 0 0 



(E-l) 

The left hand side reduces to the following form with non- 



dimensional integrand: 

4c pv(T - T ) (A 
p K v w °° 4 



1/4 



/ u*0dy* 



(E - 2) 



Since there are two distinct regions over which this integral 
is defined, it must be split into two integrals separated by 
the point, y^gg* where the inner and outer expansions should 
match. Substituting the appropriate forms for u* and 0 as 
defined by the inner and outer series the integral (E-2) can 
be written as the sum of the following two integrals: 

Inner region y* = 0 to y* = y , 

4c n v P(Tw' T J(r- ) 1/4 / Y { £ 'CnK + [f ’ Cn) g (n) + f’(n)]C 2 

P w 00 4 00 00 1 

+ [f 1 (n) g Cn) + g (n) f 1 Cn) + f’(n)]C 3 } dy* (E-3a) 
01 01 2 

where Y is the matching point, 

p = , and dy* = £dn which transforms (D-3a) into the 

following : 

_ _ r r !/4 y /F 

4c vp(T -T )(A / {f'C 2 + [f'g + f ' ] C 3 + 



w 



0 0 
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CE- 3b) 



[f'g + g £’ + f ' ] £ 3 }dn 
0 1 0 1 2 



where y depends on x. 

For the outer region, y* = y to y* = °°, the following integral 
applies : 



P i, C 

4c p pv(T w -T) (^)' 4 /°°{F ’ H + (F'H* + F’H)£ 2 

2 



r D. 



1 D 
1^2 



+ [ t ( rl ‘ 2y * } + J F ' H + i - 2y*) F"H] £ 3 }dy* (E-4) 

2 2 

The right hand integral in (E-l) also must be divided into 
two distinct regions, the first extending from the leading 
edge to the discontinuity and the second from the discontinuity 
to any arbitrary location above it. The following two inte- 
grals are obained: 



/Lq"dx + /^q” dx 



o P 



C E - 5 ) 



^Gr^ r l 

' W oo' 



where q^ is the heat flux in the undisturbed free convection 
region of the plate and q p is the heat flux above the discon- 
tinuity, both evaluated at the surface of the plate. Upon 
making suitable substitutions to non-dimensionalize the in- 
tegrands, the following result is obtained from (D-4) : 

1 

- ” -h 
^ ” 0 " 

+ r g’C 5 ] (E-6) 

2 

where Gr is evaluated at L and the g£ are evaluated at the 
surface of the plate, n = 0. 

After completing the indicated integrations in (E-3a) and 
(E-3b) and equating coefficients of like powers of £ the 
following relationships are obtained: 



- «k(T„-T JC^r [f H 1 (0) * l *ig; E‘ 
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H' (0) - H* (y ) = 3PrF(y)H(y) 



(E-7) 



C C 2 F ' (y ) H(y ) = [F"(0)] 2 y 2 



(E-8) 




D 



F'(y)H(y) + j y F ' (y ) H (y ) = - i^H'CO) (E-9) 



2 



12a 



2 



The relationships can be identically satisfied only at the 
surface of the plate, n = 0. From this fact it can be concluded 
that the outer profiles are the better approximation of global 
behavior and will satisfy the energy balance in an integrated 
sense, pointwise local behavior is required near the surface 
of the plate the inner expansions should be used, although 
theycannot satisfy the integral energy balance. 
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